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The  cross  sections  for  the  elastic  scattering  of  electrons  from 
atomic  oxygen  are  calculated  in  the  energy  range  between  0  and  13.6 
electron  volts.   Polarization  of  the  atom  is  included  by  using  a  term 
in  the  potential  function  which  has  the  correct  asymptotic  behavior  and 
which  vanishes  at  the  origin.   Previous  theoretical  and  experimental 
results  are  discussed. 
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1.   Introduction. 

A  knowledge  of  the  elastic  scattering  cross  sections  for  electron- 
oxygen  collisions  is  of  practical  importance  in  the  study  of  the  pro- 
perties of  air  at  high  temperatures.   The  experimental  determination  of 
low-energy  cross  sections  is  especially  difficult  in  the  case  of  atomic 
oxygen  because  it  recombines  easily  to  form  molecular  oxygen.   Only  re- 
cently have  experiments  been  done  to  measure  this  cross  section  in  the 
energy  range  below  15  electron  volts.  [11,15]   Several  previous  theo- 
retical treatments  have  been  made.   [1,3,7,9,12,16]   The  results  are 
divergent,  and  none  agrees  satisfactorily  with  the  experiments. 

The  elastic  scattering  of  low-energy  electrons  by  neutral  oxygen 
is  treated  here  using  the  partial-wave  method  of  Faxen  and  Holtzmark.[ 13] 
This  method,  which  is  outlined  in  the  next  section,  allows  one  to  de- 
termine the  cross  section  if  the  potential  field  in  which  the  electron 
finds  itself  is  known.   The  first  approximation  to  the  scattering  treats 
the  atom  as  a  positive  charge  at  r  =  0  surrounded  by  a  spherically  sym- 
metric negative  charge  distribution.   The  charge  density  of  this  elec- 
tronic distribution  is  calculated  from  the  Hartree  wave  functions  for 
oxygen . 

This  first  approximation  is  inadequate  because  it  ignores  the  dis- 
tortion of  the  atomic  electron  distribution  by  the  Coulomb  field  of  the 
incoming  electron.   This  polarization  effect  is  included  in  this  calcu- 
lation by  representing  it  with  a  polarization  potential  in  addition  to 
the  undistorted  potential.   The  polarization  potential  can  be  thought  of 
as  the  term  which  represents  the  difference  between  the  true  interaction 


potential  and  the  undistorted  potential.   The  form  it  assumes  for  large 
separations  is  known;  for  small  separations  the  form  is  uncertain  and 
will  be  discussed  further  in  Section  3. 

When  the  incident  electron  is  within  the  electron  cloud  of  the 
atom  there  is  an  additional  repulsive  force  due  to  the  Pauli  principle. 
Hammer  ling,  Shine,  and  Kivel  have  treated  this  exchange  effect  by  the 
inclusion  of  a  third  potential  to  approximate  the  change  in  the  electron's 
energy  due  to  the  deficiency  of  electrons  with  like  spin  near  it.  [7] 
Such  an  approximation  is  similar  to  Slater's  treatment  of  exchange  in 
atomic  wave  function  calculations.   The  result  of  their  cross  section 
calculation  which  includes  this  approximation  shows  that  the  exchange 
effect  makes  a  negligible  contribution  to  the  cross  section  in  the  energy 
range  considered  here.   Furthermore  our  polarization  potential  will  have 
a  parameter  which  is  adjusted  to  fit  experimental  data.   The  value  of 
this  parameter  is  affected  by  the  exchange  contribution,  so  no  specific 
account  is  taken  of  the  exchange  potential  in  this  calculation. 

The  actual  evaluation  of  the  cross  sections  using  the  phase  shift 
method  involves  using  a  computer  to  integrate  Schroedinger 's  equation 
numerically.   The  details  of  the  procedure  are  given  in  Section  4. 

The  following  atomic  units  are  used: 

unit  of  mass  =  m  =  rest  mass  of  the  electron 

unit  of  charge  =  e  =  charge  of  the  proton  in  Gaussian  units 

2    2   2 
unit  of  length  =  a   =  h  /4tt  me  =radius  of  the  first  Bohr  orbit 

2 
Consistent  with  these,  the  unit  of  energy  is  e  /a  ,  which  is  twice  the 

ionization  energy  of  hydrogen  and  the  unit  of  angular  momentum  is  h/2rr. 

It  should  be  noted  that  one  atomic  energy  unit  equals  about  27.2  electron 

volts;  this  is  twice  the  value  of  another  commonly  used  unit,  the  rydberg, 

which  is  13.6  ev. 


2.   Review  of  the  partial  wave  method  for  computing  the  cross  section. 

The  scattering  of  two  bodies  when  considered  in  the  center  of  mass 
system  is  equivalent  to  the  scattering  of  a  single  particle  with  the  re- 
duced mass  from  a  potential,  V(r).   For  the  electron-atom  problem  the 
reduced  mass  is  very  nearly  equal  to  the  mass  of  the  electron.   In  this 
case  the  center  of  mass  frame  is  identical  to  the  lab  frame. 

Since  the  potential  function,  V(r),  is  assumed  to  be  spherically 
symmetric,  the  electron's  wave  function  may  be  separated  and  written  as 

Y(r,0)  -2  RjCr)  P^  (cos  0). 
In  atomic  units  R.(r),  the  radial  wave  function,  satisfies  the  equation 


A_  -1 

2   dr 


^-  J  +  [2E  -  2V(r)  -  4i*±U  ]  R^(r)  «  o. 


If  we  assume  that  there  is  a  distance,  r  ■  a,  beyond  which  V(r)  is  zero, 
equation  (1)  may  be  solved  exactly  in  the  region  r  >  a,  giving 

R^r)  =  A^   [cos  ^j^Ckr)  -  sin  6£n^(kr)  ],  (2) 

where  A.   and  6.  are  constants  of  integration,  j.  and  n.  are  the  spherical 


Bessel  and  Neumann  functions,  respectively,  and  k  =V2E.   The  asymptotic 
form  of  equation  (2)  is 

>  A2 

Vr)  ~     kr~     [sin   (kr   "    *  lTJ  +  V    ]' 

so    that 

>  H 

Y    (r,9)  ~  2  P£(cos   0)  ~-  [sin   (kr  -   \  Xtt  +  b J    ].  (3) 

Ju 


•7- 


The  physical  meaning  of  6.  may  be  seen  by  realizing  that  if  there  had 
been  no  potential  anywhere,  the  asymptotic  form  of  R. (r)  would  be 

— —  sin  (kr  -  \   4tt)  .   Hence  6.,  the  phase  shift  angle,  is  the  quantity 
which  describes  the  scattering. 

In  a  symmetrical  scattering  situation  we  expect  an  asymptotic  wave 
function  to  be  in  the  form  of  a  plane  wave  superimposed  on  a  spherical 
scattered  wave.   If  the  plane  wave  propagation  is  in  the  z  direction,  we 
express  the  asymptotic  form  of  the  wave  function  as 

Y~eikz  +1  f(9)    eikr      . 
r 

i  k-z 
Since   the  asymptotic    form  of  e  is 

00 

ikz   „    .4  ...      s  _  t            s    sin  (kr-%  4rr) 
e    ~  S   i   (24+1)  P.  (cos  9)  *r * , 

4=0         *  kr 

the  asymptotic  scattered  wave  function  may  be  written  as 

Y  ~  2   i\u   +1)  P,(cos  0)  8ln  <kr^  *n)  +  f(e)eikr  .       (4) 

1=0  *  kr 

Equating   (3)    and    (4),   writing   the   sine   function  in  exponential    form, 
and   solving  for   f(9),   we  get 

00 

f(9)    =tt-    E        (24   +  l)(e2i6i-l)    P,(cos   9)    . 
Zlk     4=0  * 

The  differential  scattering  cross  section  is 

f=|f(9)|2. 
Integrating  this  over  all  solid  angles,  we  get  the  total  cross  section 

00 

a  -f§-  E   (24  +  1)  sinV  .  (5) 

k    4=0  l 
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The  value  of  6.  will  be  very  small  for  i  >   ka,  so  that  for  the  range 

of  energies  considered  here  values  of  &   greater  than  2  or  3  need  not  be 

considered  in  calculating  the  cross  section  from  equation  5. 

Hence  in  order  to  find  a,  we  need  the  values  of  6   .   These  are 

i 

found  by  solving  the  wave  equation  (1)  numerically  to  evaluate  R«(r) 
and  its  derivative  at  r  =  a.   R*(r)  and  its  derivative  are  evaluated 
at  r  =  a  from  equation  (2).   These  two  wave  functions  are  matched  at 
the  distance  a  by  equating  the  logarithmic  derivatives.   This  gives 
the  values  of  the  phase  shift  angles,  6..   When  all  significant  phase 
shift  angles  have  been  calculated,  equation  (5)  is  used  to  evaluate  the 
cross  section  at  the  given  value  of  energy.   This  process  is  repeated 
over  the  energy  range  of  interest  to  find  the  spectrum  of  cross  sections 
versus  energy. 


3   Potential  function  for  electron- oxygen  scattering. 

A  calculation  of  the  electron- oxygen  elastic  scattering  cross 
section  is  routine,  although  tedious,  if  V(r)  is  known  and  a  computer 
is  available.   In  this  calculation  we  assume  V(r)  is  the  sum  of  the 
undistorted  atomic  potential,  V„(r),  and  the  "polarization"  potential, 

n 

V  (r),  which  attempts  to  represent  the  difference  between  the  true  po- 
tential and  the  undistorted  potential. 

V„(r)  may  be  found  in  a  direct  electrostatic  manner  if  the  electron 
n. 

charge  distribution  is  known.  [14]  The  charge  distribution  is  calculated 
from  the  Hartree  wave  functions  of  the  oxygen  electrons.   Since  the 
Hartree  approximation  treats  each  atomic  electron  as  moving  in  a  spheri- 
cally symmetric  potential,  the  atomic  electron  wave  functions  can  be 
written  as 

Unim=Rn/r')Yem(9'-0'>' 

where  Y.  (0',0')  are  the  normalized  spherical  harmonics.   The  values  of 
j&m 

R  .  (r1)  have  been  tabulated  by  Hartree.  [8]   In  computing  the  electron 

* 

volume  charge  density  one  sums  UU   for  each  electron.   For  filled  n,jfc 

shells  this  sum  will  be  spherically  symmetric;  the  only  asymmetry  will 
be  contributed  by  the  outer  shell  which  is  not  filled.   It  is  within  the 
Hartree  approximation  to  neglect  this  asymmetry.   Hence  the  charge  density 
of  electrons  is  given  by 

^'>=S   NnXRL(r'> 

n,i 

where  N  „  represents  the  number  of  electrons  in  the  n,Z   shell.   Only  the 
nX 

3 
P  ground  state  of  oxygen  is  considered  here,  for  which  N   =  2,  N   =  2, 


10    '   20 


and  N   =4, 
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V  (r)  consists  of  a  negative  term  due  to  the  nucleus  and  a  posi- 
rl 

tive  term  due  to  the  electron  cloud.   It  may  be  written  as 

=°  2rr  tt 

„(,)..  i+f  [    f   [2R0(r')+2R0(r')+4R^(r')1        .,..,.„., 

H        r   J   J  J  ~ — r'  sin  0'  dO'dfl'dr' 

0  0  0         r  -  r' 


Letting  r'R^Cr')  =  ^(r'), 


"  [p2  (r«)  +  p2  (r')+  2P2  (r1)  1 
V„  -  -  -^ ^-Tn dr  . 

The  integration  can  be  divided  into  the  contribution  of  the  electrons 
between  r'  =  0  and  r'  =  r,  and  the  contribution  of  those  from  r'  =  r 
to  r'  =».   Then  V  may  be  written  as 

rl 

8       r  [P10(r')+  P20(r')+  2P21(r,)  ] 

Vu(r)  =  -  *  +   8tt  f  —& ^ ^ dr'  (6) 

H        r       J  r 

J0 


+8n  r[pio(r,>+p?>,)+2p^r'>  i 


/ 


2 
The  function  S   N  .  P  .(r1)   is  called  the  radial  charge  density.   It 
.   n*   nX 
n,Z 

represents  the  probability  per  unit  length  of  finding  an  electron  in  the 
spherical  shell  between  r'  and  r1  +  dr ' .   Figure  1  is  a  graph  of  this 
function  for  oxygen. 

The  treatment  of  the  polarization  potential  is  less  exact.   The  pro- 
cess of  treating  the  atomic  potential  as  the  superposition  of  an  undis- 
torted  potential  and  a  polarized  dipole  potential  is  less  than  rigorous. 
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However,  even  if  we  accept  this  approach,  the  form  of  the  polarization 
potential,  V_(r),  is  not  certain. 

If  the  incoming  electron  is  at  large  r,  its  electrostatic  field  at 

2 
the  atom  is  nearly  homogeneous  and  has  a  value  of  1/r  .   This  field  in- 
duces a  dipole  moment  which  is  directly  proportional  to  the  field, 

2 
M<  =  a /r  ,  where  a   is  the  polarizability .   This  dipole  moment  then  exerts 

3       5 
an  attractive  force  back  on  the  electron  with  a  value  of  f  =  2p«/r   =  2a/r 


So,  for  large  r, 


Vn(r)  -  -  f   2a/r5  dr  = 


This  form  of  V  (r)  for  large  r  must  be  modified  so  that  the  function 
does  not  blow  up  at  the  origin.   At  the  origin  there  should  be  no  polari- 
zation force  since,  in  this  limit,  spherical  symmetry  prevails  and  no 
multipole  moments  can  be  induced.   For  this  reason  the  proper  polari- 
zation potential  should  have  the  property  that 

Limit    _d    (    m 
r  -  0   dr  Pvw 

When  the  scattering  electron  comes  "inside"  the  atom,  the  atomic  elec- 
trons still  tend  to  hide  behind  the  nucleus.  However  in  this  case  the 
asymmetry  gives  rise  to  a  repulsive  force.  In  this  classical  picture, 
the  polarization  potential  should  also  have  zero  slope  at  a  value  of  r 
which  is  characteristic  of  the  size  of  the  atom. 

A  function  which  satisfies  these  conditions  and  which  will  be  used 
in  this  calculation  is 

Vp(r)  -  -  •£$     (  1  -  exp  [-  (r/rQ)q  ]  ]  .  (7) 
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This  function,  with  q  equal  either  five  or  eight,  was  first  used  by 

Biermann  to  account  for  polarization  in  the  calculation  of  oscillator 

line  strengths.  [2]   Garrett  and  Mann  have  used  this  polarization  potential, 

with  q  ■  8   in  calculations   of  electron  scattering   from  alkali  atoms.    [5,6] 

They  obtained  good  agreement  with  experimental  scattering  data  by  treating 

a   and  r  as  parameters.   Graphs  of  equation  7  are  given  in  Figures  3  and  4. 

The  exponential  in  equation  (7)  may  be  expanded  so  that  for, small 

r  we  have 

—  V  (r)  oc   r^"5 
dr  ?K   J   r-0 

From  this  we  see  that  the  requirement  that  the  polarization  force  vanish 

at  the  origin  restricts  q  to  be  greater  than  5.   In  this  calculation  we 

use  values  of  5,  6,  and  8  for  q.   The  result  with  q  =  5  represents  a 

limit  since  any  value  of  q  larger  than  five  satisfies  the  condition 

mentioned  above. 

The  value  of  r  is  taken  to  be  the  expectation  value  of  the  radius 

of  the  valance  electron  shell.   Its  value  is  calculated  from  the  Hartree 

wave  function  for  the  outer  shell  as  follows 

00 


:  =  /  r«  P^(r')  dr'. 
°   Jo    21 


This  yields  r  =1.2  atomic  units, 
o 

No  direct  experimental  values  of  the  polarizability  are  available. 
Several  calculations  based  on  indirect  experiments  have  been  made.   Klein 

and  Brueckner  found  a   to  be  5.589  by  adjusting  it  to  give  the  correct 
binding  energy  for  0  .  [9]  Using  a  similar  procedure,  Cooper  and  Martin 
calculate  a  value  of  5.499.  [3]  Both  these  methods  do  not  attempt  to 
provide  the  true  polarizability  of  the  atom,  but  yield  a  number  which 


-13- 


includes  exchange  and  other  effects.   Dalgarno  and  Parkinson  found  a 

value  of  6.005  based  on  experimental  determinations  of  the  dielectric 

constant.  [41   In  our  treatment  r   is  fixed  and  a   is  taken  to  be  either 
J  o 

5.5  or  6.0,  and  for  both  of  these  values  the  cross  section  is  computed 
for  q  =  5,  6,  and  8. 
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4.   Numerical  Methods. 

Before  the  wave  equation  can  be  integrated  it  is  necessary  to  know 

V  (r)  .   This  function  was  calculated  from  equation  (6)  using  a  trapizodial 
H 

summation  technique.   It  was  found  that  the  resulting  function  could  be 

fitted   closely  by 

1  2  -A.   r  -A   r 

VH(r)    =  -    8  IVf-^—     +Y2^~  +  Y3e  +  Y4  e  ]  , 


where  v      =   1.2806  y      =  -0.2806 

k     =   2.2376        X2  =  18.263 

Y3  =  3.0715        y4  =  -1.9710 

X-  =  6.803         X,  =  11.548. 
3  4 

This  functional  form  was  much  easier  to  enter  into  the  computer  than 
tabulated  data  points,  and  the  fitted  function  matches  the  computed   *» 
function  the  fifth  significant  figure. 

To  solve  equation  (1)  the  substitution  Y.(r)  =  r  R.(r)  was  made, 
yielding  the  equation 


d\(r)     1(1  +  1) 

\ =  [  M*  |  l)   +   2VR(r)  +  2Vp(r)  -  2E  ]  Y  (r).         (8) 

dr  r 


Equation  (8)  was  integrated  using  a  Runge-Kutta  approximation.   In  order 
to  use  such  a  method  it  is  necessary  to  know  an  initial  value  of  both  Y 
and  dY/dr  at  some  starting  point.   These  values  were  found  by  assuming 
a  solution 

00 

Y  (r)  -  va     S    C  rn  ,  C  =    1. 
s  _   n     '  o 

n=0 
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Since  this  solution  will  be  used  at  a  small  value  of  r,  a  series  approxi- 
mation for  V  (r)  and  V  (r)  may  be  used  in  equation  (8).  Neglecting  terms 
in  r  or  smaller,  we  have 


VR(r)  <*  bQ/r  +  bx  +  b2r, 


where 


bQ  =  8  (Yl  +  y2)  =8;   bx  =  8  (-  y^  -  Y2^2  +  Y3  +  Y4)  ! 


b2  =  8 


q-4 
o/  r 
Also  for  small  r,  V  «a  — —  .   Only  for  q  =  5  will  this  make  a  con- 


tribution to  the  total  potential  in  the  above  approximation,  changing  the 

a 
value  of  b  given  above  to  b  +  —  ,   Equation  (8)  for  small  r,  becomes 

1  2r 

o 


r    T* 

dr 


-  i(X  +  1)  +  2(b  r  +  b  r2  +  b  r  +  2E  r2) 


Y  =0.  (8') 
s 


Substituting  the  series  solution  for  Y   into  equation  (8'),  we  get  the 
indical  equation  o/(a   -  1)  =  Z(Z   +  1).   The  only  acceptable  root  is 
a  =  i   +  1  since  Y(0)  =  0.   The  recursion  relation  is 


C   = 
n 


■2T  b0Cn-l+(bl+E)Cn-2+b2Cn-.  1 
n(n  +  U   +  1) 


c       „/%          A  +  l            8              4+2        ,   dYs(r)         ,f  ,_>    X.    8U+2)      A+l 
So     Ys(r)^r  -  ~l     r  and  — m  (l+l)r  +  Sff*  t 
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One  could,  of  course,  evaluate  these  equations  at  r  =  0  and  start  the 

integration  from  there,  however  it  was  found  that  the  accuracy  of  the 

Runge-Kutta  integration  approximation  is  very  sensitive  to  the  step 

size  chosen  near  the  origin.   For  this  reason  we  evaluate  Y  and  dY  /dr 

s       s 

at  r  =  .01  and  use  these  values  to  start  the  iteration  procedure.   The 
step  size  was  determined  by  choosing  an  arbitrary  value  to  do  a  trial 
calculation,  reducing  the  step  size  and  repeating  the  calculation  until 
there  was  no  significant  difference  between  two  successive  results.   A 
step  size  of  .002  was  used.   Although  this  value  is  smaller  than  needed 
for  large  r,  no  attempt  was  made  to  change  the  step  size  as  the  iteration 
progressed . 

Equation  (8)  was  integrated  to  the  final  point  r  =  a  which  was 
determined  by  the  condition 

V  (a)  +  V  (a)       _4 

— — ^ =  10 

E 

For  r  >  a  we  have  a  region  of  almost  zero  potential,  so  from  equation 
(2)  the  solution  of  the  wave  equation  is 

Vr)  =  V  t  cos  bih^kr)  "  sin  6xn4(kr)  1  • 

The  derivative  is 

dY>(r) 

: =  A„r 


Vr) 


ir  2iVr  /[Xjx-i(kr)  -  w  +  u  Wkr)  i cos  6« 

-[in^Ckr)  -  (i   +  l)ni+1(kr)  ]sin  6   \  +  — 
Letting  Y^(a)  and  Y^  (a)  denote  the  computed  values  of  Y,  (r)  and 
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dr 


at  r  =  a,  we  equate  logarithmetic  derivatives  at  r  =  a  to  get 


Y^'(a)  k 

T =  

Y^(a)     (22+l)[cos  6^  ^(ka^sin  6 


^(ka)]  ( 

-[jenx_1(ka)-(X+l)nj6+1(ka)]sin6il4 


Solving  for  6.  gives 

'(ufr  *W">  -  w+d^+iO")  )  -777  Vka)]  +  Vka) 

V  Y„(a)  / 


6 .    =  arctan 


tf'Ca) 


Y,(a) 


a(2^   ^.l(ka)    "    Ci+Dn1+l(k.)    ]   -J—  n   <ka)]    +  n    (ka) 
C  Y,(a)  ) 


From  equation  (9),  for  a  given  I    and  E,  a  principle  value  of  6.  is  found, 
For  any  E,  this  process  is  repeated  for  I   =0,  1,  2,  and  3.   Then  the 
value  of  a      is  found  from  equation  (5) .  As  expected  the  contribution  to 
the  cross  section  from  the  Z   =   3  term  is  negligible  in  the  energy  range 
considered. 

The  entire  calculation  was  done  for  energies  between  0  and  .50  in 
steps  of  .05  for  q  =5,  6,  and  8;  and  for  &  =   0,  5.5  and  6.0.   All  com- 
putations were  done  on  a  C.D.C.  1604  computer  using  Fortran-60.   The  com- 
putation time  for  each  run  (a  given  q  and  a)   was  about  four  hours.   This 
large  amount  of  computer  time  could  be  reduced  by  increasing  the  step 
size  at  large  values  of  r. 


(9) 
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5.   Previous  calculations  and  experimental  data. 

The  theoretical  treatments  of  the  low-energy  scattering  problem 
may  be  divided  into  two  groups.   One  approach  is  similar  to  the  one 
used  here,  but   different  forms  of  the  polarization  potential  are  used. 
The  other  treatment  is  much  more  elaborate.   It  attempts  to  solve  the 
Hartree-Fock  equations  of  the  system  of  scattered  electron  plus  atom. 
Figure  7  shows   some    results  of  previous  calculations. 

Robinson  uses  no  polarization  potential.  [12]  His  potential  function 
is  calculated  from  Hartree-Fock  wave  functions,  and  the  cross  sections 
calculated  by  the  partial-wave  method  for  I   =  0  and  1.   He  also  evaluated 
the  cross  sections  using  the  W.K.B.  approximation  and  showed  that  it  leads 
to  very  large  cross  sections  at  low  energies. 

Klein  and  Brueckner  also  treat  the  problem  using  the  partial-wave 
method.  [9]   They  account  for  polarization  with  the  potential 

V„(r)  = 


P         2    2  2 
2(r  +  x   V 

o 

This  potential  function  has  been  used  often  in  the  treatment  of  polari- 
zation in  atomic  wave  function  calculations.   It  gives  a  zero  polarization 
force  at  the  origin,  but  differs  from  the  function  used  here  (equation  7) 
in  that  it  produces  an  attractive  force  for  all  r.   As  mentioned  before, 
they  determined  the  polarizability  from  experimental  data  on  the  binding 
energy  of  0  .   Their  curve  in  Figure  7  includes   s  and  p  waves  (i  =  0  and 
1) .   Cooper  and  Martin  [3]  found  a  mistake  in  the  paper  by  Klein  and 
Brueckner.   They  use  the  same  procedure  and  arrive  at  cross  sections  which 
show  a  decrease  at  low  energies . 
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Hammerling,  Shine,  and  Kivel  treat  polarization  in  the  same  way  as 
Klein  and  Brueckner.   They  take  o/  =  5.4  and  r  as  an  adjustable  para- 
meter. [7]   They  also  include  a  term  in  the  total  potential  to  represent 

3/2 
the  exchange  energy.   This  term  varies  as  1/r  '    and  is  derived  from 

the  Slater  approximation.   They  show  results  obtained  by  neglecting  both 
polarization  and  exchange,  and  neglecting  polarization  only.   The  cross 
sections  obtained  by  using  just  the  Hartree  potential  are  different  from 
any  obtained  previously.   The  cross  sections  with  exchange  and  polari- 
zation are  plotted  in  Figure  7 .   It  is  not  clear  why  these  should  be  so 
different  from  the  other  calculations  since  they  reported  using  the  same 
method. 

Temkin  includes  polarization  effects  within  the  Hartree-Fock  for- 
malism. [16]  He  associates  with  each  atomic  electron's  wave  function 
a  polarized  part  which  depends  on  the  coordinate  of  the  scattered  elec- 
tron.  The  wave  function  of  the  scattered  electron  is  made  antisymmetric 
with  respect  to  the  total  atomic  wave  function.   This  gives  rise  to  the 
kind  of  polarization  potential  considered  here,  called  direct  polari- 
zation, and  another  term,  called  exchange  polarization,  which  contributes 
about  15%  to  the  cross  section.   His  results  for  s-wave  scattering  are 
shown  in  Figure  7.   Bauer  and  Browne  have  modified  Temkin 's  method.  [ l] 
They  made  more  drastic  simplifications  in  the  Hartree-Fock  equations  and 
are  able  to  calculate  s,  p,  and  d-wave  scattering.   Their  cross  sections 
are  small  at  low  energies  and  fairly  flat  at  higher  energies. 

The  most  accurate  experimental  data  is  that  of  Neynaber,  et  al.  [ll] 
They  have  measured  the  cross  sections  in  the  energy  range  from  about  .1 
to  .45  (2.3  to  11.6  eV) .   They  find  the  cross  section  to  be  almost  constant 
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at  19.5  +  1.3  throughout  the  energy  range.   This  represents  a  band  in  which 
half  of  their  64  experimental  points  fall.   Sunshine,  Aubrey,  and  Bederson 
measured  the  cross  sections  between  about  .03  and  .45.   [15]  There  is 
more  scatter  in  their  data,  but  there  is  general  agreement  between  the 
two  experiments.   Lin  and  Kivel  obtained  a  low  energy  point  at  about  .02 
energy  units  using  a  shock  tube.  [10]   If  this  point  is  accurate  it  in- 
dicates that  the  cross  section  decreases  sharply  at  low  energies.   Al- 
though these  experiments  are  extremely  difficult  and  have  low  precision, 
they  do  indicate  the  magnitude  of  the  cross  section  and  show  its  general 
dependence  on  energy . 
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6.   Results. 

The  numerical  results  of  our  calculation  are  given  in  Tables  1,  2, 
and  3  and  summarized  in  Figures  5  and  6.   The  curve  labeled  a   =  0  repre- 
sents the  cross  sections  obtained  by  ignoring  the  polarization  potential. 
This  curve  corresponds  almost  exactly  to  that  of  Robinson.   This  indicates 
that  the  differences  between  Hartree  -  Fock  wave  functions  used  by  him 
and  Hartree  wave  functions  used  here  have  an  insignificant  effect  on  the 
cross  section. 

The  introduction  of  the  polarization  potential  increases  the  phase 
shift  angles  but  decreases  the  cross  section.   The  smaller  values  of  q 
produce  potentials  which  give  rise  to  smaller  cross  sections,  and  also 
straighten  the  curve  out.   The  remarkable  feature  indicated  by  Neynaber's 
experimental  results  is  the  almost  constant  value  of  the  cross  section 
between  .1  and  .45  energy  units,  [ll]  Our  results  give  a  similar  constant 
value   in  this  energy  range.   Further,  the  measurement  of  Lin  and  Kivel 
gives  evidence  that  the  cross  section  has  a  very  small  value  at  the 
lowest  energies.  [ 10]  This  is  also  a  property  of  our  results. 

Of  all  the  theoretical  calculations  that  have  been  discussed,  that 
of  Bauer  and  Brown  is  the  most  pure  (not  dependent  of  experimental  para- 
meters) which  extends  of  the  entire  energy  range,  [l]  Our  calculation 
with  <y   =  6  and  q  =  6  almost  duplicates  their  results.   Our  results 
are  also  consistent  with  those  of  Temkin  at  the  lowest  energies.  [16] 
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7.   Conclusions  and  acknowledgements. 

No  attempt  has  been  made  to  justify  the  fundamental  validity  of 
the  form  of  the  polarization  potential  except  in  a  general  way.   Experi- 
ments indicate  that  consideration  of  polarization  in  any  scattering  cal- 
culation is  important.   Although  the  "correct"  way  to  treat  the  effect 
has  been  outlined  by  Temkin,  it  is  desirable  to  have  a  more  tractable 
semi-phenomenological  method.   Our  results  indicate  that  the  potential 
function  used  here  (equations  6  and  7)  with  <x   =  6  and  q  =  5  gives  cross 
sections  which  have  the  energy  dependence  indicated  by  experiment  and 
by  the  more  exact  calculations.   Our  results  seem  much  more  realistic 
than  previous  treatments  using  other  forms  of  the  polarization  potential, 

The  author  would  like  to  thank  Professor  William  B.  Zeleny  for 
suggesting  this  subject  as  a  thesis  topic.   The  help  provided  by  the 
staff  of  the  Computer  Facility  of  the  Naval  Postgraduate  School  is 
greatly  appreciated. 
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TABLE  1 


ALPHA 


5.5 


DELTAO 


DELTA1 


DELTA2 


DELTA3 


SIGMA 


.0050 
.0500 
.1000 
.1500 
.2000 
.2500 
.3000 
.3500 
.4000 
.4500 
.5000 


9770 
8872 
8269 
7813 
7440 
7122 
6843 
6594 
6371 
6167 


1.5979 


.5018 
.5104 
.5338 
.5561 
.5759 
.5935 
.6091 
.6230 
.6356 
.6471 
.6575 


.0000 
.0055 
.0117 
.0188 
.0267 
.0354 
.0448 
.0550 
.0658 
.0771 
.0888 


.0000 
.0017 
.0035 
.0053 
.0072 
.0091 
.0112 
.0132 
.0154 
.0177 
.0200 


6 
15 
19 
21 
22 
23 
24 
24 


6678 
7451 
4121 
5178 
7739 
5334 
0109 
3357 
24.5858 
24.8064 
25.0224 


DELTAO 


DELTA1 


DELTA2 


DELFA3 


SIGMA 


.0050 
.0500 
.1000 
.1500 
.2000 
.2500 
.3000 
.3500 
.4000 
.4500 
.5000 


9709 
8719 
8078 
7598 
7209 
6878 
6591 
6335 
6105 
5897 
5706 


.5006 
.5274 
.5611 
.5889 
.6121 
.6318 
.6488 
.6638 
.6771 
.6890 
.6998 


.0000 
.0055 
.0117 
.0188 
.0266 
.0353 
.0447 
.0547 
.0653 
.0765 
.0880 


.0000 
.0017 
.0035 
.0053 
.0072 
.0091 
.0112 
.0133 
.0154 
.0177 
.0201 


10.4912 
22.2733 
27.5908 
30.0143 
30.9885 
32.2318 
31.1230 
30.8598 
30.5507 
30.2516 
29.9895 


DELTAO 


0  -    8 
DELTA1 


DELTA2 


DELTA3 


SIGMA 


.0050 
.0500 
.1000 
.1500 
.2000 
.2500 
.3000 
.3500 
.4000 
.4500 
.5000 


9644 
8554 
7871 
7366 
6960 
6617 
6320 
6057 
5822 
5610 
5416 


.5028 
.5582 
.6009 
.6318 
.6562 
.6762 
.6934 
.7083 
.7214 
.7332 
.7438 


ALL  PHASE  SHITT  ANQbES  AP6  QIV6N 
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.0000 
.0055 
.0117 
.0187 
.0265 
.0351 
.0444 
.0543 
.0648 
.0757 
.0869 


in  MybTiPbes  er 


.0000 

15 

,9588 

.0017 

36, 

,8699 

.0035 

42 

,9553 

.0053 

43 

,8549 

.0072 

43, 

,0616 

.0091 

41, 

.7041 

.0112 

40, 

,2216 

.0133 

38, 

,7979 

.0154 

37, 

,5062 

.0177 

36, 

,3703 

.0201 

35, 

3891 

er  pi 

RADIANS 

TABLE  2 


ALPHA  =  6.0 


.0050 
.0500 
.1000 
.1500 
.2000 
.2500 
.3000 
.3500 
.4000 
.4500 
.5000 


DELI A0 


9825 
8977 
8388 
7939 
7569 
7253 
6976 
6729 
6505 
6301 
6114 


0  =  5 

DELTA1 

.5023 
.5040 
.5234 
.5434 
.5619 
.5786 
.5937 
.6074 
.6198 
.6312 
.6416 


DELTA2 

.0000 
.0060 
.0129 
.0207 
.0294 
.0390 
.0496 
.0609 
.0729 
.0856 
.0987 


DELTA3 


SIGMA 


.0000 

4, 

,0100 

.0019 

12, 

,8560 

.0038 

16, 

,3795 

.0058 

18, 

,5432 

.0079 

20, 

,0096 

.0100 

21, 

,0657 

.0122 

21, 

,8826 

.0145 

22, 

,5658 

.0184 

23, 

,1803 

.0193 

23 

,7635 

.0218 

24, 

,3341 

DELTAO 


.0050 
.0500 
.1000 
.1500 
.2000 
.2500 
.3000 
.3500 
.4000 
.4500 
.5000 


9756 
8810 
8181 
7707 
7320 
6992 
9670 
1.6450 
1.6220 
1.6012 
1.5820 


DELTA1 

.5014 
.5186 
.5483 
.5745 
.5970 
.6165 
.6336 
.6487 
.6622 
.6744 
.6854 


DELTA2 

.0000 
.  0060 
.0128 
.0206 
.0293 
.0389 
.0493 
.0605 
.0724 
.0848 
.0976 


DELTA3 

.0000 
.0019 
.0038 
.0058 
.0079 
.0100 
.0122 
.0145 
.0168 
.0193 
.0219 


SIGMA 

• 

7.4262 
18.2930 
23.2639 
25.9655 
27.4278 
28.2089 
28.6273 
28.8643 
29.0231 
29.1604 
29.3036 


DELTAO 


.0050 
.0500 
.1000 
.1500 
.2000 
.2500 
.3000 
.3500 
.4000 
.4500 
.5000 


9685 
8632 
7959 
7458 
7054 
6712 
6415 
6152 
1.5917 
1.5703 
1.5509 


DELTA1 

.5006 
.5426 
.5830 
.6140 
.6390 
.6599 
.6778 
.6934 
.7072 
.7196 
.7307 


DELTA2 

.0000 
.0060 
.0128 
.0206 
.0292 
.0387 
.0490 
.0600 
.0717 
.0838 
.0963 


DELTA3 


SIGMA 


.0000 

12, 

,2868 

.0019 

28, 

,7581 

.0038 

35, 

,5747 

.0058 

37, 

,9080 

.0079 

38, 

,3704 

.0100 

38, 

,0360 

.0122 

37, 

,3904 

.0145 

36, 

,6634 

.0168 

35, 

,9643 

.0193 

35, 

,3414 

.0219 

34, 

.8102 
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TABLE  3 


ALPHA  =  0.0 


DELT-AO 


DELTA1 


DELTA2 


DELTA3 


SIGMA 


.0050 
.0500 
.1000 
.1500 
.2000 
.2500 
.3000 
.3500 
.4000 
.4500 
.5000 


1.9154 
1.7502 
6672 
6121 
5705 
5370 
5091 
4850 
4638 
4450 
4280 


.5004 
.5157 
.5493 
.6000 
.6664 
.7415 
.8139 
.8750 
.9223 
.9575 
.9836 


ALL  PHASE  SHIFT  ANGLES  ARE  GIVEN  IN 


.0000 

.0000 

86.6421 

.0000 

.0000 

63.6804 

.0004 

.0000 

51.5073 

.0011 

.0000 

48.9037 

.0020 

.0000 

53.3931 

.0032 

.0002 

60.4892 

.0047 

.0004 

64.6529 

.0064 

.0005 

63.9188 

.0084 

.0008 

59.9304 

.0105 

.0011 

54.7727 

.0128 

.0014 

49.6356 

MULTIPLES 

1  OF  PI 

RADIANS 
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